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ABSTRACT 

Using the observations of the star formation rate and H I densities to z ~ 4, with measurements of 
the Alolecular Gas Depletion Rate (MGDR) and local density of H2 at z = 0, we derive the history 
of the gas consmiiption by star formation to z 4. We find that closed-box models in which H2 
is not replenished by H I require improbably large increases in p(i?2) and a decrease in the MGDR 
with lookback time that is inconsistent with observations. Allowing the H2 used in star formation 
to be replenished by H I does not alleviate the problem because observations show that there is very 
little evolution of phi{z) from z = to z = 4. We show that to be consistent with observational 
constraints, star formation on cosmic timescales must be fueled by intergalactic ionized gas, which 
may come from either accretion of gas through cold (but ionized) flows, or from ionized gas associated 
with accretion of dark matter halos. We constrain the rate at which the extraglactic ionized gas must 
be converted into H I and ultimately into H2. The ionized gas inflow rate roughly traces the SFRD: 
about 1-2 xlO* M© Gyr~^ Mpc~'^ from z ~ 1 — 4, decreasing by about an order of magnitude from 
z = 1 to z = with details depending largely on AIGDR(t). All models considered require the volume 
averaged density of PH2 to increase by a factor of 1.5 - 10 to z ^ 1.5 over the currently measured value. 
Because the molecular gas must reside in galaxies, it implies that galaxies at high z must, on average, 
be more molecule rich than they are at the present epoch, which is consistent with observations. These 
quantitative results, derived solely from observations, agree well with cosmological simulations. 
Subject headings: Galaxies:ISM — Galaxies:evolution — Stars:Formation 



1. INTRODUCTION 

The time variation of the mean star forma- 
tion rate in galaxies is by now well established 
dMadau ct al. 1996 ; Steidel et al. IQJI iHippelein et al.l 
120031 : iHopkins fc Beacomll2006D . The star formation rate 
is either flat or slowly rising with time, reaching a maxi- 
mum near z ~ 1 - 2, and then declines precipitously down 
to the current epoch. This change in the star formation 
rate must be closely coupled to both the inventory of gas 
available for star formation, and the way in which this 
gas is channeled into galaxies. 

Stars condense only from molecular gas at the cur- 
rent epoch and at all epochs in the past. This state- 
ment derives from both observational and theoretical 
considerations. Observationally, the youngest stars are 
always found to be associated with their nascent molec- 
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most completely molecular, representi ng a stable phase o f 
the ISM with little atomic content (iBurton et al.lll978D . 
Theoretically, there is general consensus that the initi- 
ation of star formation requires the nascent gas to be- 
come Jeans unstable, probably mediated by magnetic 
fields (|Shu et al.l I1987D . In star forming regions, T is 
typically 10 - 20 K, but in any event cannot be less 
than 2.7 K. In order to get a solar mass star at 10 K, 
the Jeans instabiliy criterion would require a density of 
Pj > (kT / ^imuG fiTT^ l^&M]) ~ lO^cm^^ if a molecular 
core forms a star at 100% efficiency. The density must be 
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higher if the efficiency is lower as many observations now 
sugge st (e.g. iMotte et al.]ll998l : lAlves et al.l 120071 : iMversI 
120081 ). In order to reach these temperatures and densi- 
ties, the gas must be fully molecular in order to achieve 
the necessary cooling. We would therefore expect that, 
in a broad sense, the gas consumption rate is closely tied 
to the star formation rate. 

Unfortunately, there are few constraints on the gas 
from observations because little is known about the dis- 
tribution of neutral gas at high z. There are very few 
detections of atomic gas in emission at z > 0.1, and 
what little we know about the atomic gas comes from 
Lyman-alpha lines seen in absorption towar d quasars 
and radio-loud AGN (e.g. IP^o chas'ka & Wolfe 200l 
iWolfe et all [20051 : iZwaan fc Prochaska .2006) . Molecu- 
lar line observations at high-z have been largely lim- 
ited to the brightest objects, though some recent obser- 
vations at z ^ 2 have begun to probe lower luminos- 
ity systems (iFo rster Sc hreiber et all 120091 : iDaddi et al.l 
[200alTacconiet al.. 20101 1 

Gas consumption by star formation in galaxies has 
been investigated previously via observations of the gas 
depletion time, Tdep = Mgas/ SFR. This represents the 
amount of time it will take the galaxy to completely ex- 
haust its gas supply at the current star formation rate. 
Studies of individual local disk galaxies find depletion 
times on the order of a few Gyr, much shorter than the 
Hubb le time (e.g. iLarson et al.l Il980t iKennicutt et al.l 
Il994f ). This is the gas depletion problem: without some 
form of gas replenishment, star formation in disk galaxies 
should be coming to an abrupt end. One propo s ed so - 
lution is stellar recycling, which iKennicutt et al.l ()1994l ) 
find can extend the gas depletion times in many local 
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disk galaxies by a factor of 1.5 to 4. Gas accretion has 
long been suggested as a solution to the gas depletion 
problem as wel l, and observational support for inflow is 
accumulating. iSancisi et al.l ()2008[ ) review the observa- 
tional evidence for gas accretion such as galaxy inter- 
actions and minor mergers, high velocity clouds (HVCs), 
extra-planar gas and warped or lop-sided HI disks. These 
observations yield an estimate for the 'visible' gas accre- 
tion rate onto a typical disk galaxy of ~ 0.2 Mq yr^^, 
which falls short of the typical star formation rate of ^ 1 
M0 yr~^. These studies focus on recycling and inflow in 
local disk galaxies, but one may ask how does gas con- 
sumption evolve with redshift? 

The time evolution of gas in disk galaxies has been 
studied on large scales using observations of damped 
Lya systems (DLAs) to infer the cosmic mass density 
of HI as a function of redshift. The cosmic mass den- 
sity is the mass density averaged over a large volume so 
as to be representative of a typ ical Mpc~^ of th e uni- 
verse. iLanzetta et"all (|1995D and lPei fc Falll (|1995D build 
simple models of gas consumption focusing on the chem- 
ical evolution of the gas using the cosmic mass density 
of HI and observed metallicities as inputs. Their models 
predict the star formation rate density (SFRD) of the 
universe as a function of redshift for their chosen inflow 
and outflow parameters. However, the SFRD is now be- 
coming increasingly well constrained by observations and 
we can use SFRD(z) as an input to our models in order 
to place constraints on the gas inflow rates. 

Finally, this problem has been approached using nu- 
merical simulations to explain observations of the SFRD 
and galaxy properties. The shape of the SFRD has been 
investigated using cosmological, hydrodynamical simula- 
tions that include prescriptions for star formati on and 
feedback (e.g. Hcrnquist & Springel 2003; Schave ~t al.l 
120101 ) . These studies suggest that the shape of the SFRD 
at high redshifts is determined by the buildup of dark 
matter halos and the gas brought in with them, and the 
decline of the SFRD at low redshifts is due to lower cool- 
ing rates in the halo gas, gas exhaustion and stellar and 
black hole feedback. The results of simulations have also 
been used to build simple models in order to better un- 
derst and gas accre t ion to fuel star formation in galax- 
ies. iBouche et al.l ()2009D build a simple model of gas 
consumption starting with simulated halo growth his- 
tories with different prescriptions for gas accretion and 
star formation. The predictions from the different pre- 
scriptions are compared to the observed SFR-Mass and 
TuUy-Fisher relations for galaxies from z ^ 2 to z = 0. 
The authors find that the prescription for gas accretion 
modeled on cold flows with halo mass cutoffs agrees best 
with observations. 

In this paper, we look at the issue in reverse. We build 
simple models of gas consumption based solely on ob- 
servations in order to understand the roles of the differ- 
ent phases of gas in star formation in galaxies. Using 
observations at 2; = 0, we examine which quantitative 
conclusions can be extrapolated to high 2, and using ob- 
servations of the star formation history, we make sev- 
eral inferences about how the gas consumption into stars 
must have proceeded. We find that the relationship of 
the inventory of gas to that of the stars is not straight- 
forward: the observations imply that that all phases of 
the interstellar and intergalactic medium must be taken 



into account in order to understand how the gas forms 
stars. Building a model that includes all the gas phases, 
we make predictions about gas densities and their varia- 
tion at intermediate- and high-z, and how this gas must 
have been accreted into galaxies. 

The paper is organized as follows. 321 describes the 
observations of the SFRD, MGDR, and pHi that we 
use as inputs to our models. In fj3l we build three models 
to fit the observations: the restricted closed box model, 
the general closed box model and the open box model. 
In fJH we discuss potential changes to our star formation 
prescription at high redshift and examine the predictions 
of the open box model. Throughout this paper, we adopt 
a standard ACDM cosmology with (h, Qm, ^a) — (0.7, 
0.3, 0.7). All of the densities are in comoving units. 

2. OBSERVATIONS 
2.1. SFRD 

The observed (comoving) star formation rate den- 
sity, SFRD, was estimated by iMadau et all (|1996[ ) over 
a large range of z and re-examined by several other 
investigators (iSteidel et al.lll999t fHippelein et "all 120031: 
Hopkins & Beacom"2006l) . We will focus on the results of 
Hopkins & Bcacom (2006(), a compilation of SFRD mea- 
surements at different wavelengths. The measurements 
are converted to a common cosmology, SFRD calibra- 
tion, and dust obscuration correction; the data are then 
fit to a piecewise linear form in log(l-l-z) vs. \ or(SFRD) 
space as well as the parameterizat ion from ICole et all 
(|200in . IHopkins fc BeacomI (|2006f ) find that changing 
the assumed IMF corresponds to a simple change in the 
overall amplitude of the SFRD, so each fit is done us- 
ing two extreme IMF forms in order to provide bounds 
on the actual form. These two extreme forms for the 
IMF are a modified Salpcter A IMF and the form of 
IBaldrv fc G lazcbrook (2003). Fig.djplots the fits for the 
two IMFs in red and blue. We have smoothed th e origi - 
nal piecewise linear fits from IHopkins fc BeacomI (|2006D 
for our models (solid lines). The form of the smoothed 
piecewise linear fit is given in Appe ndix A. 

At z = 0, the SFRD fits from IHopkins fc BeacomI 
(pool predict a local SFRD of (0.8 - 1.5) x 10"^ 
Mq Mpc~^ yr~^. In an extens ive study of galax ies in 
the local universe (to z = 0.1), ISalim et al.l mOlh find 
that the SFRD is 1.828l°;J;^^ x 10" ^ hro Mq Mpc- ^ yr"! 
from UV observations (using the IChabrieil 120031 IMF) , 
a value they argue is the most accurate determination 
of this number to date. This value (shown as a data 
point in Fig. 1) is inde ed close to those found from other 
recent investigations fHouc k et al.l 120071 : iHanish et al.l 
2006) and is also in reasonable agreement with the forms 
of iHopkins fc BeacomI |2006). 

The agreement between the locally determined 
SFR D and extrapolat io n fro m studies at higher 
z (IHopkins fc Beacon] 120061: ISteidel et afl 119991: 
iHippelein et al.l I2003D " is encouraging and provides 
some confidence that the star formation rates and 
their time variation are being measured reliably. While 
somewhat different functional forms for the decline 
of the SFRD with time have been proposed in the 
literat ure, e.g., log ( SFRp ) has been found to be linea r 
in z (jSteidel et al]|1999D, in t CHippclcin et al.l l2003l) . 
and in log(l -|- z) (|Hopkins fc Bcacom .20()6() . these 
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Fig. 1. — Comparison of forms of the SFRD. The sohd 
curv es show our smoothed version of the piecewise linear fits 
from IHo pkin s fc Beacora |2003) ; the dashed curves show the fits 
to the Cole ct aT] 1120011 ) form. Two IMFs are plotted; the 
modified Salpeter A IMF (SalA IMF; red) and the IMF from 
IBaldrv GlazebrooH l|2003l) (EG IMF; blue). The green s ymbol 
denotes the SFRD at z = measured bv lSalim et al.l 1)20071 '). 

differences have only a small quantitative effect in what 
follows. 

2.2. MGDR 
2.2.1. Measurements at z — 

The star formation efficiency, SFE, is often defined as 
the star formation rate per comoving volume divided by 
the mass per com oving volume of gas; its units are yr~^ 
(jLerov et al.ll2008l and references therein) . Defined in this 
way, the SFE is not properly an efficiency, but a rate, 
and we drop this unfortunate usage even though it has 
become firmly entrenched in the observational literature. 
Since we are interested in SFEm, the star formation rate 
density divided by the density of molecular gas, , we 
introduce the molecular gas depletion rate, MGDR, to 
replace the usage of SFEm- We define MGDR as SFRD 
divided by the density of molecular hydrogen, (this 
does not include He). The inverse of MGDR is the molec- 
ular gas depletion time, tm, which represents the time it 
takes to consume all of the molecular gas at the current 
rate of star formation. 

Recently, iLerov et al.] ()2008[ ) have measured MGDR(z 
= 0) from a comprehensive analysis of 23 nearby galax- 
ies. The results are based on H I su rface densities mea - 
sured from the THINGS H I survey (IWalter et al.l[2Q08l) . 
H2 surface densit ies inferred from the BIMA SONG 
(|Helfer et al.ll2003n and HERACLES (Leroy et al. 2008b) 
CO surveys, and star f ormation rates f rom both SINGS 
jKennicutteLalJ |2003l) and GALEX (|Gil de Paz et al.l 
|2007) data. The galaxies surveyed include spiral and 
dwarf galaxies and the analysis was done on a pixel- 
by-pixel basis convolved to a common resolution, typi- 
cally '^SOO pc. This comparison is the most extensive 
work done to date and the authors find a remarkable 
constancy of MGDR over a wide range of conditions: 
0.525 ± 0.25 Gyr~^, equivalent to a molecular gas de- 
pletion time of 1.9 Gyr. Their measured star formation 
rates and H2 column densities vary by three orders of 
magnitude averaged over entire galaxies, and the pixel- 
by-pixel values vary even more. Thus the constancy of 



the MGDR occurs over a wide range of conditions both 
within galaxies (including nuclei and disks) and from 
galaxy to galaxy. 

2.2.2. Measurements at high z 

The ILerov et"ari (|2008f ) study only applies to galaxies 
near z = 0. To extend the range of z we appeal to 
observations of the MGDR in normal, z > 1 galaxies. 
[Daddi ct al. ( 2009) report the SFR and total gas mass for 
six, near-infrared selected BzK galaxies at z ~ 1.5. Using 
numerical simulations, they calculate a conversion factor 
aco = Mgas/Lco = 3.6 ± 0.8 M© (K km s^^ pc^)"!. 
This value is close to a Milky Way-like value of aco "^4.6 
and excludes a typical ULIRG value of aco ^ 1- This 
gas mass includes He, so we divide by 1.4 to calculate 
the MGDR, which we have defined to not include He. 
The M GDR values for these six galaxies in IDaddi et al.l 
(I200I v ary between 1.9 an d 4.7 Gyr ^, three to nine 
times the lLerov et al.l (j2008[) value for z = 0. Apparently, 
molecular gas is consumed by star formation at a much 
higher rate at high redshift than it is today. 

Similarlv. lTacconi et al.l ()2010f ) have made an extensive 
survey of normal, star forming galaxies at redshifts 1 and 
2, measuring the MGDR for a sample of 19 galaxies. At 
each redshift locus, the selected galaxies sample the high 
mass end of the main sequence galaxy pop ulation in the 
M*-SFR plane. The iTacconi et all (|201O) resuhs use a 
Milky Way-like value for aco- The y find MGDR values 
in agreement with the IDaddi et al.l (f2009f ) values at z 
1.5. The data points from these st udies are plotted in 
Fig ures Blan d [5l Da ddi et al.l ()2009[ ) as purple diamonds 
and lTacconi et al.l (j2010[r as green squares. 

In samples of gal axies selected by d ifferent optical 
and near- IR criteria, iReddv et all ()2005[) find that BzK, 
BX/BM and DRG selected galaxies account for an SFRD 
of - 0.1 Mq Mpc"^ yr^i in the range 1.4 < 2 < 2.6. This 
is most of t he observed SFRD (see Fig. [T]). More recently, 
i Reddy e"t~a l. (2008) find that galaxies with Ljr w Lboi ^ 
10^2 Lq account for « 70% of the SFRD at 1.9 < z < 2.7. 
The BX/BM, BzK and DRG selection criteria sample 
luminous, star forming galaxies w ith L ^ 10^^ — 10^ ^ 
Lq and SFR ~ 10 - 500 M© yr"! fTacc oni et al1l2008l) 
The refore, these systems sampled by Daddi et al.l ()2009l) 
and ITacconi et al.l (|2010f ) account for most of the SFRD 
at these redshifts so that the MGDR value that typifies 
these galaxies should describe the average cosmic MGDR 
in Eq. [1] This motivates us to use these samples to con- 
strain our guessed forms of the MGDR at other redshifts. 

Numerous other authors have estimated the MGDR 
from individual hi gh redshift galaxies, or from starbursts 
or ULIRGS (e.g. iGao fc Sololm^ I200I . and all have 
found that the MGDR is greater in these galaxies than 
is typical of galaxies at 2 = 0. In fact, there is no obser- 
vational evidence for a declining MGDR with increasing 
redshift to at least z '--^ 4. Taken together, all of these 
observations lead us to reject any model of gas consump- 
tion that requires lower MGDR values at redshifts signif- 
icantly greater than zero. 

2.3. 

Using a combina tion of CO and H I measureme nts in 
the local universe, lObreschkow &: Rawlingsl ([2009') have 
determined the density of H2 at 2 = 0, Ph2{0)j to be 1.9 - 
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2.8 xlO^ h Mq Mpc~^. The range of values corresponds 
to different assumptions about how metallicity affects the 
CO-to-H2 conversion ratio (the authors quote an error of 
±40% on each individual calculated value). Using h — 
0.7 gives (0) = (1.3-2.0) x 10^ Mq Mpc'^, or an aver- 
age value of 1.65 x 10^ Mq Mpc~^. This value is within 
50% of the estimate of Ph.(O ) = 1-1 x lO"^ MgMpc-^ 
by iZwaan fc Prochaskal (|2006l ) ( no error quoted) using 
a different set of observations and should therefore be 
reasonably reliable. 
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points fo r th e comoving pjfj from 

Rao et al.l pOpg i: ILah et al.l 0UTi ): and 

Proch aska &: Wolfe ( 2009 ) with hnear fit in Iog(l + z) vs log(pjj/) 
space. 



2.4. PHI 

A number of recent studies have reported measure- 
ments of the mean H I comoving mass density in galax- 
ies, pht{z). from t he local universe to z ^ 6. At z = 0, 
IZwaan et al.l ()2005[) present a study of the H I mass func- 
tion using the H I Parkes AU Sky Survey (HIPASS) data. 
At low redshift, Lah et al. (2007) calculate pni by co- 
adding H I 21-cm emission from galaxies with known po- 
sitions and redshifts. At higher redshifts, es timates of 
Pm are mainly obtained from DLA studies (|Rao et al.l 
[200l iProchaska fc Wolfi [20091 ) . We use these observa- 
tions of H I to construct an analytical form for pni as 
a function of z. We fit a straight line to these points in 
log(l -I- z) vs \og{pHi) space. The data points and the 
fit are shown in Fig. O The observations suggest very 
little evolution of pui , which increases by only a factor 
of 2-3 between z = and z = 6. A more recent DLA 
study (|Noterdaeme et al.l 120091 ) using SDSS DR7 finds 
the values at z ^ 2 t o 3 to be somewhat higher than the 
IProchaska fc Wol^ (|2009| ) points; this would make our 
linear fit in Fig. [2] correspond even better to the obser- 
vations. 

3. BUILDING A MODEL TO FIT THE OBSERVATIONS 

The observations above allow us to compare the evo- 
lution of the SFRD and dpn^/dt. We start with the sim- 
plest model possible: a closed box of H2 being turned into 
stars at the rate of dp/f^/dt. In this closed box model, 
we initially assume that the MGDR is constant in time 



(the restricted closed box model; ij3.ip and subsequently 
relax this condition (the general closed box model; ii3.2p . 
The failure of both models motivates us to consider an 
open box model ( §3.3|) . where we allow the densities of 
all four phases of the IGM - stars, H2, H I and H II - to 
vary to match the observational constraints. 

We assume that the molecular gas is depleted only 
through star formation, and that any H2 dissociated or 
ionized by star formation is instantaneously returned to 
the molecular state. Given the short timescales for the 
formation of molecular gas from its atomic form, y 10^ 
yr at the relevant densit ies (jHollenbach fc SalpetedlT97ll : 
iCazaux fc Tielenj|2004D , this approximation should be a 
good one for the purposes of this paper. In any event, 
if we define dpH2 / dt to be the net fiow rate of molecular 
gas into stars, then there is no ambiguity. 

We write the statement that star formation occurs only 
through the depletion of molecular gas: 



SFRD = MGDR x pn^ 



(1) 



This equation was used to infer the individual MGDR 
for a sample of nearby galaxies with a wide of range 
of SFRs and H2 column densities, where it was found 
that MGDR « 0.5 Gyr~^ to a remarkable constancy 
(|Lerov et al.ll2008L see §2.2|) . We use this prescription for 
the SFRD for our models but we note that the form may 
change at higher redshift (See t j4.ip . 

It is also interesting to note that if we divide the ob- 
served global star formation rate den sity SFRD(z = 0) ~ 
(0.8 - 1.8) X 10-^Mq M pc"'^ yr-i (jHopkins fc BeacomI 
[200I iSalim et al.|[20Ci7l) by the observed pg,(0) = (1 . 3 - 
2.0) X IO^Mq Mpc-3 ()Obreschkow fc Rawling!j[2009l ) . we 
obtain a range of MGDR that is consistent with the val- 
ues of iLerov et al.l (|2008[ ). Given that stars must form 
from molecular gas, this result is not surprising. Never- 
theless, the agreement is reassuring because it is based 
on different data sets and different methods of determin- 
ing the relevant quantities. It also suggests, combined 
with the arguments in fJT]that Eq. [T]can be extrapolated 
to all z. 

We use mass densities instead of mass surface densi- 
ties, which are used by observers, but note that these are 
roughly equivalent because for the most part, the stars, 
H I and H2 are generally confined to thin disks within 
galaxies. 

3.1. The Restricted Closed Box Model 

In the closed box model, we consider only stars and 
II2, and allow pH2 to be converted into stars at the star 
formation rate density, SFRD: 



dpH2 
dt 



-SFRD. 



(2) 



For the moment, we consider a restricted closed box 
model in which we take the MGDR to be constant as a 
function of redshift. To assess the ability of this model 
to fit observations, we first combine eqs. (1) and (2) to 
obtain d{SFRD) /dt = MGDR x SFRD. We then note 
that our piecewise linear fit to the observed SFRD as a 
function of time (Appendix A) implies d{SFRD)/dt ^ 
(.24Gyr^^) x SFRD, where the coefficient is about half 
of the observed MGDR (,Lerov et al...200a ). 
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In other words, from the assumption of a constant 
MGDR at the present epoch in a closed box H2 model, we 
find that the star formation rate is declining only half as 
fast as expected given our current reservoir of molecular 
gas. It is this discrepancy in the derivatives of the ob- 
served cosmic star formation rate and the rate at which 
we observe molecular gas being converted into stars that 
we call the cosmic molecular gas depletion problem. We 
note here that given the uncertainties in the observa- 
tions, this factor of two in itself is not a strong argument 
against the closed box model, but we will show in gen- 
eral that observational constraints rule out any closed 
box model. 

3.2. The General Closed Box Model 

We now relax the assumption of a constant MGDR in 
the closed box model in tj3.1l To study the predictions 
of this model, we calculate (z) by integrating Eq. ([2]) 
and using the observed SFRD(z) as an input. We then 
divide SFRD(^) by ph2{z) to obtain MGDR(z). The 
results are shown in Fig. Bl where ph-, (0) is set to the 
mean value from lObreschkow fc Rawlingsl (j2009() . The 
uncertainties in the SFRD (due to the IMF) discussed in 
^2.1\ are seen to have only a minor effect on the resulting 
PH2iz) and MGDR(z). 

Fig. [3] shows that pH2 (t) increases by a factor of ^ 10 
from z = to 1, and MGDR decreases with increasing 
redshift, contrary to the observational results discussed 
in ii2.2\ Thus, even the general closed box model is at 
odds with the observations, leading us to our next model. 

3.3. The Open Box Model 

Since a closed box model of only II2 and stars is incon- 
sistent with observations, we now allow additional com- 
ponents that can be converted into H2 and then into 
stars. We consider separately the H I gas and an exter- 
nal source of gas that we call pext, and modify Eq. ([2|) 
to 

^^^^ ^SFRD (3) 
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dt 



dt 



dt 



3.3.1. 



Phi 



For the H I gas, the observations discussed in §2.41 and 
Fig. [5] suggest that pHi{z) is very slowly varying over 
cosmic timescales; therefore dpni/dt is small. Fig. 2] 
shows that the derivative of the observed pni (red curve) 
is an order of magnitude smaller than the observed SFRD 
(black curves). In the absence of Pe^t, we have \dpH2/dt\ 
(blue curves) is approximately equal to SFRD, as in the 
failed closed box model. Thus the inclusion of H I alone 
in an open-box model is not enough to fit the model to 
the observations. This leads to a robust conclusion: the 
reservoirs of HI and H2 at all times in the past (at least 
as far back as z = A) are insufficient to fuel the star 
formation over cosmic timescales. 

It is important at this point to clarify that 
Phi represents the reservoir of H I both in galaxies as 
well as the H I outside of galaxies. This is because the 
DLA observations of Prochaska and Wolfe (2009) include 
all of the high column density neutral H I (N(HI) > 2x 
10^° 2). This gas contains at least 85% of the neutral H I 
atoms for z < 6 (O'Meara et al 2007). 
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Fig. 3. — The predicted MGDR(t) and pH2 (*) in the closed box 
modeL The top panel shows the input SFRD forms (the smoothed 
piecewise linear fits from Fig. [!}• The middle and bottom panels 
show the MGDR and pjj^ (t), respectively, predicted by the general 
closed box model for these observed SFRD(i) forms. The points in 
the m iddle panel show the MGDR o bservations fro m [Daddi et al] 
12009) (purple diamo nds) and from ITacconi et al.l 12010 ) (green 
squares), discussed in il2.2l This model requires a lower MGDR in 
the past, contrary to observations. 



3.3.2. peait 

We are therefore forced to include a nonzero dp^xt/dt 
term in the open box model. This component represents 
the ionized intergalactic gas at all temperatures and den- 
sities that can recombine to form H I within a Hubble 
time. Effectively, it is the ionized gas in the filaments of 
the cosmic web. 

Note that some of the H I can become ionized and re- 
distributed to the intergalactic medium at high enough 
temperatures and low densities such that this gas does 
not recombine in a Hubble time. Such gas can be ejected 
by means of supernovae, galactic winds or AGN. How- 
ever, if we define dpext/dt to be the net flow out of the 
ionized phase, then the gas that is fed back into the ion- 
ized state is implicitly included in our accounting. That 
is, any H I fed back into the ionized phase is made up by 
an equivalent increase in dpext/dt. If the total amount of 
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Fig. 4. — Rates of gas flowing from one phase to another in 
an open box model with only stars, H2 and HI (ignoring pext )■ 
The SFRD (black curves) is an order of magnitude higher than 
dpHi/dt (red curve), forcing dpH2/<H (blue curves) to be approx- 
imately equal to the SFRD, as i n th e closed box model. The two 
forms of the SFRD described in ^2TT1 are used: modified Sal A IMF 
(solid) and the Baldrv & Glazcbrook (2003) IMF (dashed). The 
derivatives of p^j and p^f, are negative, but the absolute values 
are plotted here. 



ionized gas available is represented by ftbaryon minus the 
total amount of baryons in galaxies, we may consider the 
ionized phase to be a nearly infinite reservoir of gas avail- 
able to fuel star formation. Any gas ionized and added 
to that reservoir by star-formation and active galaxies is 
negligible. 

To compute the dpext/ dt{t) required in the open box 
model to match the observations in §2, we start with an 
observed SFRD(t) and pHi{t), and a guess for the form 
of the MGDR(t) that is compatible with the data points 
fromlTacconi ct al. (2010) and Daddi ct al. (2009). From 
these inputs, we compute (t) and its time derivative 
using pH^it) = SFRD(i)/MGDR(t). Combining dp//, /dt 
with the observed SFRD(t) and dpHi{t)/dt in Eq. §1 
then gives dpext/dt. The results of this procedure are 
shown in Fig. [5l 

To illustrate the expected range of possible values, 
we use two forms for the input MGDR (top left panel) 
and two forms for the input SFRD (black dashed curves 
in the right panel). We bracket the possible values 
for the MGDR on one side as constant using the mea- 
sured value at z = 0, and on the other as one that in- 
creases linearl y to l arger redshifts as sugg ested by the 
iTacconi et all (|2010( ) and lDaddi eFall (|2009l ) data points. 
For the SFRD, we use the smoothed piecewise fits from 
iHopkins &: Beacoml ()2006D for the two extreme IMFs dis- 
cussed in ^2.11 We then calculate dpj/^/dt, dpext/dt and 
PH2 for each of the four possible combinations of MGDR 
and SFRD. The range of calculated values is indicated 
in Fig. [5] by plotting the minimum and maximum of the 
four curves for each calculated quantity. 

In the bottom left panel, we note that the minimum 
PH2 curve lies below pui at all times, and the maximum 
curve becomes larger than pni at z ~ 0.35. The differ- 
ence between the two curves is mostly the MGDR choice: 
the maximum curve corresponds to the higher SFRD and 
the flat MGDR, while the minimum curve corresponds to 
the lower SFR and the increasing MGDR (as expected 
from Eq. [T]). All other combinations of MGDR and 
SFRD whould lie between these two curves. Because all 



of the molecular gas resides in galaxies, we require that 
the molecular gas mass in galaxies is larger, on average, 
at z = 1-2, than is typical at z = 0. The change in pH2 
resulting from changing the SFRD is minor. Therefore, 
observations of pH2 at z > or the redshift at which 
PH2 — Phi will constrain the form of the MGDR. 

In the right panel, the absolute value of the rates is 
plotted, with a thicker line style used to indicate negative 
rates. The negative portions of the dpH2/dt curves corre- 
spond to decreasing pH2 as we move forward in time to- 
wards z = 0, as H2 is converted into stars. The dpext/dt 
curves are negative for the whole range of redshifts plot- 
ted, indicating a flow of external gas into the H I reser- 
voir of galaxies. We note that the range of solutions for 
dpext/dt does not deviate very much from the SFRD: 
roughly a factor of two at low redshifts for the minimum 
case. This is because the reservoirs of H2 and H I are 
so small compared to what is required by the observed 
SFRD. Therefore, we conclude that the amount of inflow 
needed from this external gas, dp^xt I dt, is approximately 
equal to the SFRD. This observational conclusion is rein- 
forced by cosmological simulations which find that star 
formation rates cl osely f ollow gas infal l rates (Keres et aLl 
[20051 iDekel eran[200l . Specificallv. (Keres et all (|2009l ) 
calculate an upper limit on the external gas supply feed- 
ing galaxies that is only a factor of 2 higher than our pre- 
dicted dpext /dt curves and shows similar evolution with 
redshift. 

4. DISCUSSION 



We have shown in §3. H and §3.21 whv the closed II2 box 
model doesn't work. Here we discuss the open box model 
of §3.31 and what this model predicts. 

4.1. Variations in the SFRD 

For the models in ^ we have extended the star for- 
mation rate prescription from lLerov et al.l (|2008[ ) to mass 
densities averaged over Mpc scales: SFRD oc pjy,. Some 
studies of galaxies with higher gas surface densities have 
foun d evidence for a steeper power law of SFRD oc 
p]it (|KennicuttJ [19981 : IWong fc Blit!^ [20021 : IBouche et all 
|2007| ). The choice of SFRD prescription at a given red- 
shift depends on what type of galaxies dominate the 
SFRD at that redshift. At the presen t epoch, re gular spi- 
ral galaxies, such as those studi ed bvlLerov et al . (2008), 
seem to dominate the SFRD: iSalim et al.l (|2007l) finds 



that galaxies in the mass range 10^-^ < M* < lO^^-^M© 
account for about half of the total SFRD. At higher red- 
shifts, however, the BzK and BX/BM galaxies dominate 
the SFRD (see ij2.2.2p . These galaxies are more gas rich, 
and thus the prescription may be more appropriate. 

For this paper, we have assumed the iLerov et"aLl 
(|2008| ) prescription for the SFRD, but note that the 
power law may change at higher redshift because the 
dominant mode of star formation may change. This 
is equivalent to a change in the MGDR with time, a 
case we consider explicitly and pr obably contrib u tes to 
the values of the MGDR found by ITacconi et all (pOlOl ) 
and "DaddTeTaD (120091 ). If, for example, SFRD = 
MGDR{Q)p]f^ at higher redshift, the n we would find 
MGDR{z) = MGDR{Q)p^j:il by forcing our prescrip- 
tion of SFRD = MGDR x p/f, . Thus as p/f, increases 
at higher redshift, an increase in the power law of the 
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Fig. 5. — Predictions of the open box model, calculated from combinations of the inputs: two forms of the MGDR (solid black lines in 
top left panel) , two forms of the SFRD (dashed black lines in right panel) and the fit to the observations of pu i (red dotted line in bottom 
left panel). In the top left panel, we plot the MGDR points from n^cconi et~al.l (|2010 ) (green squares) and IDaddi et aLl 1)20091 ) (purple 
diamonds) to constrain the guessed range of MGDR forms. The possible range of each of the calculated quantities, dpjj^/dt, dpext/dt and 
shown by two lines illustrating the minimum and maximum curves from the four combinations of inputs. In the right panel, the 
absolute values of the rates are plotted: thicker lines indicate negative values. 



SFRD prescription will be manifested as an increase in 
the MGDR. An increase would tend to bring the MGDR 
closer to the upper bound in the top panel of Fig. [5l with 
the result that {z) would lie closer to the lower blue 
curve in the bottom panel in Fig. [5l 

4.2. Stellar Recycling 

Stellar recycling is an important co mponent of any 
treatment of gas evolution in galaxies. iKennicutt et al.1 
(UqM), for example, suggest that gas returned to the ISM 
during stellar evolution can signficantly increase the gas 
depletion time in the Milky Way and in nearby galaxies. 
In the model presented here, we do not explicitly con- 
sider the effect of gas return to the ISM in Eqs. [T] and 
[3l but we argue that our formulation already accounts 
for stellar recycling due to the nature of the observed 
quantities we use as inputs. Most of the ret urn of gas t o 
the ISM comes from red giant stars (see e.g. lBlitdll997l ). 
and much of that is returned to the ISM in the form of 
molecules (Marengo: 2009). Even gas that is returned in 
other phases becomes largely molecular after each spiral 
arm passage, the timescale for which is ~ 10^ yr in most 
galaxies, a timescle short compared to the rather long 
timescales we consider in this paper. Since the recycled 
gas is largely molecular, and any that isn't is quickly 
returned to the molecular phase, we argue that observa- 
tions of molecular gas already include the gas returned 
via stellar recycling. Therefore, this recycled component 
is included in our initial condition, pHziz = 0), as well 
as the MGDR. Since we use the MGDR to relate SFRD 
and at each time step, stellar recycling is built into 
our model through these observations, so we need not 
include an explicit recycling term in our equations. 

4.3. Behavior of PH2 (z) 



Although the exact shape of pH2 depends sensitively 
on the form of MGDR, we have bounded the behavior 
of PH2 by calculating what we take to be limiting cases 
in our open box model (Fig. [5]). The pH2 curves all rise 
with increasing redshift, peaking around z ^ 1 — 1.5 at 
1.5 to 10 times the value of ph2 today. After the peak, 
PH2 may fall off toward higher redshift if the MGDR is 
rather constant, or stay close to constant if the MGDR 
remains high. The prediction is less we ll constrained 
at hig her redshifts. It is noteworthy that iTacconi et al.l 
([2010') find that the gas disks they observe at 2 = 1 and 
z = 2 have considerably more molecular gas relative to 
the stars, typically about 30-50%, compared to ^ 1 — 5% 
for th e Milky Way and nearby disk galaxies (jHelfer et al.l 
|2003() . This trend is consistent with our estimates. Al- 
though there is some uncertainty in the II2 masses of 
ITacconi et al.l (2010) because of the uncertainty in the 
value of Xco , there appears to be little doubt that the ra- 
tio of H2 to stellar mass in the BzK and BX/BM galaxies 
is higher than typical values for similar galaxies at z = 0. 
Due to the sensitivity of PH2 to the form of MGDR, fu- 
ture observations of PH2 or the ratio PH2/ Phi at higher 
redshift would allow us to better constrain the evolution 
of MGDR, and reduce the area between the bounding 
curves of Fig. [5l 

4.4. The Nature of dpext/dt 

In our open box model, the H I reservoir, pni, is 
augmented by an inflow of gas from p^xt at a rate of 
10^ - 10*MoMpc-3 Gyr"\ depending on z. This high 
rate of inflow means that the gas being accreted is mostly 
ionized since the fraction of neutral hydrogen outside of 
the Phi reservoir is too small. 

The Phi inferred from observations of DL A systems ac- 
counts for the H I associated with galactic disks. As men- 
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tioned in Sg^l lO'Meara eFall (|2007D find that systems 
with column density I^hi < 2 x 10^° cm~^ account for 
« 15% of neutral hydrogen atoms at all z < 6. Therefore, 
the fraction of H I outside DLA systems is about 15%, 
corresponding to roughly 1.5 x 10^ MoMpc"^. For an 
average inflow rate of a few times 10'' Mq Mpc^^ Gyr~^ 
for the past 10 Gyr, this intergalactic H I could only ac- 
count for ten percent of the total, at most. Therefore, 
the inflow of gas needed for fueling ongoing star forma- 
tion represented by dp^xt/dt must be almost completely 
ionized. 

Recently, cold flows have been suggested as an impor- 
tant source of gas for galaxy formation an d evolution 
(|Keres et al.l 120051 : iDekel fc BirnboimI l200l . In these 
models, galactic disks in halos with M < are 
built up by direct accretion of cold gaseous streams from 
the cosmic web. For galaxies with larger masses, cold 
flows are also the dominant means of mass accretion, 
but different outcomes for individual galaxies depend on 
the epoch of inflow. If this picture is correct, our work 
implies that the cold flows must be almost entirely ion- 
ized. 

The same is true if the gas needed to fuel star formation 
is brought in primarily through minor mergers. If this gas 
were in atomic form, it would be part of the inventory of 
atomic gas observed in the DLA systems, which we have 
shown in ij3.3l to contribute negligibly to fueling the star 
formation at all redshifts up to z = 4. 

4.5. Cooling Times 

The open box model requires dp^xt/dt ^ SFRD, or 
about 10^ to 10^ MqMpc"^ Gy"^ We use these num- 
bers for dpext/dt to calculate a cooling time for the gas 
in the context of two models for gas ac cretion onto galax- 
ies: t wo-phase cooling of hot halo gas (Mailer & Bulloc^ 
I2004f) and cold flow accretion (Kercs & Hernquist 200^ 1 

We estimate the cooling time, tcooi, by taking 



Pgas 
tmnl 



Pext ; 



(4) 



where pgas is the average mass density of the cooling 
ionized gas smoothed over the appropriate volume (to 
be specified for each cooling model individually). We set 
pgas equal to ruHnef where Ue is the local number den- 
sity of electrons and / is the filling factor for the relevant 
volumes (ne/rie). Combining this with the cooling time 



tc 



2A(r)ne 

where A(T) is the cooling function of the gas, gives 

ShTpext 



2A{T)mH 



(5) 



(6) 



As a basis for comparison, we first estimate the filling fac- 
tor / for hot halos of galaxies out to the virial radius. 
We make the simplistic assumption that the universe is 
made up of L^, galaxies with masses Mdyn ^ 10^^ Mq and 
circular velocities of ^ 160 km s~^. Therefore, the virial 
radius, Ryir = GMdyn/v^ ^ 300 kpc. We estimate the 
average number density in this simple universe composed 
of i» galaxies by dividing the total luminosity density, 
£, by L*. We use the r band values from .Blanton et al.i 



(|2003D which calculates the galaxy luminosity function at 
z ~ 0.1 from SDSS data: £ w 1.84 x 10^ h Mq Mpc'^, 
L, w 1.2xl0i"h-2 Mq. This yields ni. - C/L^ ~ 0.015 
h^ Mpc~^ ^ 5 X 10"^ Mpc~^. More recent work on the 
lumin osity function using SPSS DR 6 yields similar re- 
suhs (jMontero-Dorta fc Pradal |2009i). Therefore, in this 
simple universe, the filling factor for the galaxy halos 
is 

4 



/ 



nL,-TTR: 



3 

vir 



6 X 10" 



(7) 



iMaller fc Bullocfl (|2004l ) consider gas within the cool- 
ing radius of a halo, Rc, cooling via cloud fragmentation. 
This results in the formation of warm 10"* K) clouds 
within the hot gas halo. In this model for the two-phase 
cooling of the hot halo gas, the relevant temperature for 
the gas is the virial temperature of the halo, ^ 10^ K for 
a Milky Way type ha lo. 

iKeres fc Hernquist (|2009f ) find that the majority of 
cold clouds that form around Milky Way type galaxies 
are the result of filamentary " cold mode" accretion. Most 
of the gas does not reach the virial temperature of the 
halo, ^ 10^ K, but rather cools from a maximum tem- 
perature of - 10^ - 10^ K. 
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Fig. 6. — Shaded regions show the allowed filling fraction / and 
jie for three temperatures: 10^ K in blue, 10^ K in green and 10^ K 
in red. The axis on the right gives the radius of the relevant volumes 
as a fraction of the virial radius corresponding to the filling fraction 
/ in a simple universe filled with L, galaxies. At each temperature, 
the allowed region of jig vs / space is calculated using the predicted 
range for the inflow rate at z = 0: dpext/dt = 0.4 — 1.2 X 10^ 
Mq Mpc~^ Gyr~^. The region is further bounded by the vertical 
dotted line at the value of jie corresponding to a cooling time equal 
to the age of the universe. This line is only plotted for the 10 K 
gas: for the other two temperatures, the densities are very small 
and lie outside of the plotting range. 

We examine the possible values for and / in these 
two models by calculating / as a function of rig for three 
temperatures: 10"* K and 10^ K for cold flow accretion 
and 10^ K for cooling from the hot halo. In Fig. HI 
we plot / versus rig for these three temperatues for the 
range of dp^xt/dt at z = predicted by our open box 
model: 0.4 - 1.2 x lO"^ Mq Mpc'^ Gyr"!. The y-axis on 
the right side shows \og{R/ Rmriai) corresponding to /, 
where R is the radius of the relevant volume associated 
with an galaxy. The vertical dotted line indicates the 
value of rie at which the cooling time is equal to the age 
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of the universe. We use the approxim ate form for A(T) 
for mi ldly enriched gas {Z — 0.1) from lMaller fc Bullockl 
((2OOI : 



A(r) ~ 2.6 X 10^23 



T 
106 K 



-1 



cm erg s ^ 



(8) 



For the cold flow gas at 10'' K and 10^ K, the dotted 
lines where the cooling times equal the age of the universe 
are outside of the plotting range, so the minimum allowed 
value of Ue is where R ~ Rviriai- about 8 x 10"'' cm~'^ 
for the 10^ K gas and even smaller for the 10'' K gas. For 
the cooling hot halo gas at 10^ K, the condition that the 
cooling time be less than the age of the universe forces 
rie to be larger than about 2 x 10~^ cm~'^. 

These calculations don't put strong constraints on the 
density of the halo gas since the cooling times are so rapid 
for a large range of temperatures and densities. How the 
gas gets into the galaxies themselves will involve a more 
complete treatment including the effects of self-sheilding, 
which is beyond the scope of this paper. 

4.6. Comparing dpext/dt to Dark Matter Accretion Rate 

The rate dpext/dt inferred from our open box model 
provides an estimate for the average rate at which the 
baryonic fuel is required to make its way to a galac- 
tic disk in order to sustain the observed star formation, 
which largely occurs in the disk. A comparison of this 
rate with the mean rate of baryon accretion at the virial 
radius of the host dark matter halo will provide an es- 
timate for the efficiency of converting the cosmological 
infalling baryons into stars. Many cooling and feedback 
processes obviously affect the fate of baryons after their 
infall onto the halo and whether they will reach the disk. 
In fact, much of the current research in galaxy formation 
modeling is aimed at understanding this transition. Our 
goal here is to estimate an overall ratio, as a function of 
redshift, of the baryon accretion rates at the virial radius 
and at the disk scale. 

We begin with the dark matter accre- 
tion rate from iMcBride et al.l ()2009f ) and 
iFakhouri. Ma fc Bovlan-KolchinI (j2'om ). which quan- 
tified the mass accretion histories of all dark matter 
halos with masses above ^ lO'^" M0 in t he two Millen- 
nium simulations of a ACDM un iverse (jSpringel et al.l 
120051: iBovlan-Kolchin et al.l I2009D . An approximate 
function is provided for the average mass accre- 
tion rate as a function of redshift and halo mass 
(jFakhouri. Ma fc Bovlan-KolchinllMol ): 



M = /3M0yr-i 



M 



(9) 



where 7) — (25.3, 1.65) for the median rate and (46.1, 
1.11) for the mean rate, D,„i and i^A are the present- 
day density parameters in matter and the cosmological 
constant, and flm + ^^a is assumed to be unity (as in 
the Millennium simulation). The mean rate is ~ 50% 
higher than the median rate due to the long tail of halos 
with high accretion rates in the distribution. This M 
represents the average rate at which the mass in dark 
matter is being accreted through the virial radius of a 



halo. The mass growth comes in two forms in cosmolog- 
ical simulations: via m ergers with other resolved halos 
(jFakhouri fc Mall2009aD . and via "diffuse" accretion of 
non-halo material that is a combination of unresolved ha- 
los and unbound dark matter particles (jFakhouri fc Mai 
,2009b '). 

We convert M above into a mean accretion rate for 
the baryons, Mh, by assuming a cosmic baryon-to-dark 
matter ratio of = rif,/f2,„ = 1/6. The result should 
provide a reasonable approximation for the mean rate of 
baryon mass that is entering the virial radius via merg- 
ers plus accretion of intergalactic gas. These infalling 
baryons are presumably in a mixture of forms: warm- 
hot ionized hydrogen gas of 10^ to lO^K, "cold" flows of 

IC^K (still ionized) gas, and H I, H2, and stars brought 
in from merging galaxies. As discussed earlier, the ma- 
jority of these baryons must be in the form of H II gas. 

To compare Mb with the rate of external gas inflow, 
Pext ( §3.3p . needed to account for the evolution of the 
observed star formation rates, we define 



dpext/dt 



fbMM{dn/dM) 



(10) 



where M is the mass of the dark matter halo under con- 
sideration, M is calculated using Eq. and dn/dM 
is the (comoving) number density of dark matter halos 
with mass in the range of M and M + dM . The parame- 
ter a represents the fraction of accreting baryons (at the 
virial radius) that must be converted into stars in our 
open box models. 

The value of a can be estimated by combining the 
allowed range of dpext/dt from Fig. [5] with the halo 
abundance dn/dM from the Millennium simulation 
(jSpringel et al.ll2005l ). Taking fb = 1/6 and M = lO^^ 
Mq, we find the predicted a to be ~ 70 — 100% at z > 3, 
- 120 - 200% at z ~ 2, and - 30 - 90% at z = 0. 
We note here that several factors may contribute to the 
uncertainty in the alpha values. First, the distribution 
of M is broad at a given h alo mass (see, e.g.. Fig. 5 
of IFakhouri. Ma fc Bovlan-K olchin 2010), and we have 
simply used the mean value for a rough estimate of a 
here. Second, a depends on the halo mass appropri- 
ate for the population of galaxies that dominates the 
SF RD at a given redshift . Using a clustering analy- 
sis, (AdelbergereEaD (EQQI) find that BzK and BX/BM 
galaxies have halo masses of about 10^2 Mq, but the aver- 
age value appropriate to our analysis may vary. However, 
since M cx M^'^ and Mdn/dm is approximately oc M~^, 
the mass dependence is weak, so alpha changes by a max- 
imum of ^ 20% if we change the halo mass by a factor 
of 3. Third, since the dpext/dt calculated in our open 
box model roughly traces the SFRD, the alpha values, 
especially around z ^ 1 — 2, will be affected by the exact 
form of the SFRD. Finally, the large alpha value at z ~ 2 
may be reflecting a change in the fraction of baryons in 
the fllaments at that redshift. Considering all the un- 
certainties, we make the conservative suggestion that the 
open box model requires a large fraction (~ 30 — 90%/ of 
the infalling baryons at the virial radius to be turned into 
stars from z ^ — 4. This is consistent with the work of 
iDekel et al.l (j2009D . which flnds that the star formation 
rates in the typical 'star-forming galaxies' at 2; ~ 2 are 
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very close to the baryonic inflow rates from simulations. 

Another way to estimate a is to compare the various 
M directly. At z 0, the median baryon accretion rate 
from cq. (9) is Mb - 9 M© yr^^ for 2 x lO^^ j^a- 
los, and the measured star formation rate in the Milky 
Way ranges from M ^ , ^ 2 M© yr"^ to ~ 4 M© yr~^ 
(IMiller fc ScaIol[T97llDiehl et al.li2006D . These rates im- 
ply a ~ M^/Mb ~ 25 to 50% at z - 0. In addition, 
we can use the predicted pext from our open box models 
to estimate a mean conversion rate of external H II gas 
into stars per galaxy, M ~ Pext/nL,- Taking p^xt ~ 10^ 
MgMpc"^ Gyr^i from Fig.[5]and n^. ~ 5 x 10"^ Mpc'^ 
for the number density of L* galaxies today, we obtain 
a rate of ~ 2 Mq yr~^, which is consistent with the star 
formation rate in the Milky Way. 



5. SUMMARY AND CONCLUSIONS 

In this paper, we have built a quantitative model 
of gas consumption on cosmic scales based solely on 
observations. Using the observed the Star Formation 
Rate Density (SFRD), Molecular Gas Depletion Rate 
(MGDR), and volume averaged density of molecular hy- 
drogen {PH2) a-nd atomic hydrogen (phi) we have defined 
the cosmic molecular gas depletion problem and calcu- 
lated the mass flow rates and densities of the star- forming 
gas back to z ~ 4 needed to resolve it. Extrapolations 
further back in time are primarily limited by uncertain- 
ties in the SFRD. We find that: 



• There are no models of gas consumption where the 
reservoir of gas is either only H2 or both H I and 
H2 that are compatible with the observations of 
molecular gas in the galaxies at z 2 that produce 
the observed SFRD. 



• Inflowing ionized intergalactic gas must provide 
most of the gas that turns into stars to times as 
early as z ~ 4. There is so little neutral gas inflow 
at all epochs, that the neutral gas ought to be con- 
sidered more as a phase of cosmic gas flow than a 
reservoir of star forming gas. 

• The rate of mass inflow from the ionized state to 
the atomic state roughly traces the star formation 
rate density. From z ~ 1 — 4, the mass inflow 
rate is 1 - 2 xlO^ Mq Mpc"^ Gyr"!. At z < 1, 
the mass inflow rate varies linearly from about 0.5 
xlO^ Mq Mpc-3 Gyr-i at z = to about 1.5 xlO^ 
M0 Mpc"3 Gyr-^ At ah redshifts, we find the 
mass inflow rate must be a significant fraction of 
the inf ailing rate at the virial radius, in agreement 
with simulations. 

• For all models, the volume averaged density of II2 
increases from its present value by a factor of 1.5 
to 10 at z = 1 - 2 depending mostly on MGDR(t). 
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APPENDIX 



A. FUNCTIONAL FORM FOR THE SFRD 



The form for the SFRD used throughout the paper is a smoothed form of the piecewise fits from lHopkins fc BeacomI 
p006, ). The original piecewise function is defined by intercept and slope a and b for z < zi and c and d for z > zi 
(we only consider z < 4). The smoothed section is just a quadratic form in the log(l + z) vs log{SFRD) plane over a 
distance 2A in log(l + z). We give the piecewise function for the SFRD below, with x — log(l + z) and xi = log(l + zi). 



log (SFRD) 



bx 

-b{xi 

-lib- 
dx 



d+ib 



4A 



d)^ 



A 



x^ - {xi - 

{x — Xl 



f A) 



Xl 
Xl 



<a;< Xl - A 
A <x< Xl + A 



(Al) 



In the table below, we give a,b,c,d and zl for each of the two piecew ise linear fits from lHopkins fc BeacomI ()2006f ) 
(for the Modified Salpeter A IMF and the IBaldrv fc Glazebrookl ()2003l ) IMF). To produce reasonable smoothing, we 
used 2A = log(1.5). 





Mod SalA IMF 


BG 2003 IMF 


a 


-1.82 


-2.02 


b 


3.28 


3.44 


c 


-0.724 


-0.930 


d 


-0.26 


-0.26 


Zl 


1.04 


0.97 



TABLE 1 

Table of fit parameter values for the two piece wise linear fits fromT Jopkins fc BEAcdI3 (120061) : the Modified Salpeter, A 

IMF and the IBaldry fc Glazebrook (.2003') IMF. 



